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LONG-TERM  GOALS 

The  long-term  scientific  goal  of  my  research  is  to  better  understand  the  distribution  of  phytoplankton 
in  the  world's  oceans  through  remote  sensing  their  influence  on  the  optical  properties  of  the  water. 

OBJECTIVES 

Optically,  phytoplankton  reveal  their  presence  through  their  influence  on  the  inherent  optical 
properties  (lOP’s)  of  the  water.  The  main  effect  of  phytoplankton  is  to  increase  the  absorption  of  light 
by  virtue  of  the  strong  absorption  by  their  photosynthetic  (chlorophyll  a)  and  accessory  pigments.  A 
secondary  effect  is  to  increase  the  backscattering  coefficient  of  the  medium  in  a  manner  that  depends 
on  the  concentration  of  pigments.  Although  techniques  for  measuring  the  absorption  coefficient 
directly  (e.g.,  in-situ  AC9  measurements  or  in-vitro  filterpad  absorption)  are  becoming  accepted  by  the 
scientific  community,  laboratory  techniques  for  measuring  backscattering  are  tedious  and  subject  to 
error,  and  in-situ  techniques  for  backscattering  are  in  their  infancy.  In  addition,  in  most  in-situ 
measurements  the  volume  of  medium  that  is  sampled  is  small  and  may  not  be  representative  of  the 
whole  water  body,  even  in  a  homogeneous  medium.  Thus,  in  the  past,  there  has  been  considerable 
effort  devoted  toward  indirectly  inferring  these  lOP’s  by  virtue  of  their  affect  on  the  apparent  optical 
properties  (AOP’s),  e.g.,  the  diffuse  reflectance  of  the  water  (the  color  of  the  water)  or  the 
downwelling  irradiance  attenuation  coefficient.  These  AOP’s  are  perhaps  the  most  frequently 
measured  quantities  in  hydrologic  optics.  Clearly,  interpretation  of  such  observations  requires  a 
detailed  understanding  of  the  influence  of  phytoplankton  on  the  lOP’s,  and  their  link  to  the  AOP’s. 

The  lOP  AOP  link  forms  the  focus  of  the  present  research.  In  particular,  our  research  is  centered  on 
deriving  the  lOP’s  from  measurements  of  the  AOP’s.  This  is  an  example  of  the  inverse  problem  of 
radiative  transfer.  It  is  important  in  that  lOP’s  determined  from  AOP’s  are,  by  definition,  sampled  at  a 
scale  appropriate  for  radiative  transfer,  and  for  remote  sensing.  Also,  the  retrieved  lOP’s  possess  the 
attribute  that  when  combined  with  the  radiative  transfer  equation,  they  reproduce  the  measured 
AOP’s.  The  table  below  summarizes  in  matrix  form  inverse  problems  relevant  to  hydrologic  optics. 
Our  goal  is  to  fill  in  the  uncompleted  elements  of  the  matrix. 
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Matrix  of  Inverse  Problems  Relevant  to  Hydrologic  Optics 


Problem 

Water  Body 

Homogeneous 

Stratified 

Elastic  (Depth  =  oo) 

✓ 

✓ 

Elastic  (Depth  <  oo) 

✓ 

+ 

Raman  (Depth  =  oo) 

X 

X 

Raman  (Depth  <  oo) 

+ 

+ 

Eluor.  (Depth  =  oo) 

Eluor.  (Depth  <  oo) 

Inverse  solution  completed  earlier 
X:  Inverse  solution  completion  reported  here 
+:  Inverse  solution  to  be  studied  during  the  next  reporting  period 

In  earlier  research  (Gordon  and  Boynton,  1997,  1998),  we  focused  on  developing  methods  for  deriving 
the  vertical  profiles  of  the  absorption  and  backscattering  coefficients  from  measurements  of  the 
vertical  profiles  of  downwelling  irradiance  and  either  upwelling  radiance  or  upwelling  irradiance  in  the 
absence  of  inelastic  processes  (^).  However,  as  the  presence  of  inelastic  processes  can  make  a 
significant  contribution  to  the  in-water  light  field  (Stavn  and  Weidman,  1988;  Gordon  1999),  we 
extended  that  study  to  include  the  influence  of  Raman  scattering  in  the  inversion  algorithm,  i.e.,  to 
retrieve  the  lOP’s  from  the  AOP’s  when  they  are  contaminated  by  Raman  scattering  (X).  We  chose  to 
look  at  Raman  scattering  as  opposed  to  other  inelastic  processes  (e.g.,  fluorescence  of  dissolved 
organic  matter)  because  the  Raman  scattering  cross  sections  are  known  (Bartlett,  et  ah,  1998). 

APPROACH 

Our  approach  for  including  Raman  scattering  in  the  water  body  is  to  use  the  basic  algorithm  developed 
by  Gordon  and  Boynton  (1998),  to  retrieve  the  absorption  and  backscattering  coefficients  of  a 
vertically  stratified  water  body.  We  refer  to  this  algorithm  as  the  elastic  inversion  algorithm  (EIA). 
However,  as  the  AOP’s  that  are  entered  into  the  algorithm  are  corrupted  by  Raman  scattering,  we  use 
irradiances  measured  at  the  Raman  excitation  wavelength  to  estimate  the  Raman  source  function  and 
then  its  contribution  at  the  emission  wavelength  of  interest.  This  contribution  is  subtracted  from  the 
irradiances  at  the  emission  wavelength  and  the  EIA  applied  to  the  remainder.  Because  the  Raman 
contribution  to  the  irradiance  at  the  emission  wavelength  depends  on  the  emission  wavelength  lOP’s, 
this  process  must  be  carried  out  iteratively,  i.e.,  new  estimates  of  the  lOP’s  are  used  with  the  source 
function  to  develop  new  estimates  of  the  Raman  contribution,  etc. 

WORK  COMPLETED 

Using  the  approach  described  above,  we  developed  an  algorithm  for  retrieving  vertical  profiles  of  the 
absorption  and  backscattering  coefficients,  in  the  presence  of  Raman  scattering  in  a  vertically  stratified 
ocean  from  measured  vertical  profiles  of  the  AOP’s  at  both  the  emission  and  Raman  excitation 
wavelengths.  Details  of  the  Raman  inversion  algorithm  (RIA)  are  provided  in  Boynton  and  Gordon 
(1999). 


RESULTS 


As  an  example  of  the  perfoiinance  of  the  RIA  we  used  the  Morel  and  Gentili  (1993)  model  and  the 
pigment  profile  in  Figme  1  to  provide  a  set  of  lOP’s  to  generate  synthetic  data.  The  pseudodata  were 
generated  at  460,  544,  and  666  mn.  We  assumed  there  was  no  Raman  contiibution  at  460  nm,  but 
Raman  scattering  from  460  to  544  nm  was  included  in  the  544  nm  light  field,  and  Raman  scatteiing 
fr  om  544  to  666  was  mcluded  in  the  666  nm  light  field.  The  data  at  666  nm  were  fust  inverted  using 
just  the  EIA  algorithm  yielding  the  tr  iangles  in  Figiue  2  for  a(z)  and  bb(z).  Note  that  the  enor  in  a(z)  is 
as  much  as  20%  and  the  enor  in  bb(z)  reaches  a  factor  of  10.  Using  the  Raman  inversion  algorithm 
(RIA),  these  enors  are  reduced  to  1-2%  for  a(z)  and  about  20%  for  bb(z).  This  improvement  is  typical 
of  the  perfoiinance  of  the  RIA  in  cases  where  the  Raman  contribution  to  the  light  field  is  lar  ge.  It  is 
interesting  to  note  that  the  Raman  contiibution  from  544  nm  to  the  total  upwelling  madiance  (666  mn) 
at  12  m  in  Figme  2  is  over  90%.  The  algorithm  breaks  down  as  it  must  when  the  Raman  contribution 
overwhelms  the  elastic  contribution,  because  it  is  then  impossible  to  estimate  the  Raman  contribution 
with  sufficient  acciuacy.  As  in  the  EIA,  an  attractive  featiue  of  the  RIA  is  that  the  results  are  nearly 
independent  of  the  scattering  phase  fimction  assumed  for  the  medium  in  the  inversion  process. 


Depth  (meters) 


Figure  1:  The  pigment  concentration  (in  mg/m^)  used  to  test  the 
Raman  inversion  algorithm. 


Retrieved  a(z)  at  666  nm 


Retrieved  bb(z)  at  666  nm 


Figure  2:  Performance  of  the  Raman  inversion  algorithm  (RIA)  in 
retrie\'ing  a(z)  and  bb(z)  compared  with  the  elastic  inversion  algorithm 
(EIA)  for  the  pigment  profile  in  Figure  1. 


IMPACT/APPLICATIONS 

We  believe  that  these  algorithm  will  be  of  significant  utility  for  processing  existing  and  fiitiue 
experimental  iiradiance  profile  data  to  estimate  the  absorption  and  backscattering  coefficients,  and 
then  relationship  to  constituent  concentrations,  for  irse  in  ocean  color  remote  sensing  algorithms. 
Fmthermore,  the  algorithm  provides  lOP's  that  are,  by  definition,  sampled  at  a  scale  appropriate  for 
radiative  transfer,  and  therefore,  will  be  of  significant  vahre  for  exanrinirrg  qirestions  of  closrne 
conceruing  traditional  lOP-instrarments  that  sample  at  scales  of  a  few  cirbic  centimeters. 

TRANSITIONS 

We  have  prepared  versions  of  the  EIA  (most  irsefitl  in  the  gr  een  and  bhre  regions  of  the  spectnrm)  to 
nm  on  work  stations  using  Linrrx  (GNU  Fortran)  or  Windows  NT  or  95/98  (Power-Station  4  Fortran). 
The  latter  has  been  provided  to  H.  Sosik  for  testing. 

RELATED  PROJECTS 

We  are  working  with  J.  Mireller  and  J.R.V.  Zaneveld  to  compare  om-  estimates  of  absorption  and 
backscattering  profiles  with  direct  measiuements  irsing  in-sitir  instnmientation.  We  expect  this  work 
will  contr  ibirte  to  AOP-IOP  closiue.  We  are  developing  a  bi-diiectional  model  of  the  water-leaving 
radiance  (inchrding  Raman  contribirtions)  for  Case  1  waters.  This  will  be  irsed  itr  processing  MODIS 
imager-y,  and  is  fimded  by  NASA. 
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